Thermal fluctuations in superconducting nanowires 
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We re-analyze the problem of thermally activated phase slips (TAPS) which can dominate the 
behavior of sufficiently thin superconducting wires at temperatures close to Tc- With the aid of an 
effective action approach we evaluate the TAPS rate which turns out to exceed the rate found by 
D.E. McCumber and B.I. Halperin, Phys. Rev. B 1, 1054 (1970) within the TDGL analysis by the 
factor ~ (1 — T /Tc)' 1 S> 1. Additional differences in the results of these two approaches arise at 
bias currents close to the Ginzburg-Landau critical current where the TAPS rate becomes bigger. 
We also derive a simple formula for the voltage noise across the superconducting wire in terms of 
the TAPS rate. Our results can be verified in modern experiments with superconducting nanowires. 



I. INTRODUCTION 

Fluctuations are known to play an important role in 
superconducting structures with reduced dimensions. In 
the case of superconducting nanowiresi fluctuations may 
essentially determine the system behavior in a wide tem- 
perature interval causing, e.g., non- vanishing wire resis- 
tance down to T = 0. 

Over four decades ago it was realized by Little 2 - that 
sufficiently thin superconducting wires may acquire a 
non-zero resistance below the BCS critical temperature 
of the bulk material Tc due to non-trivial thermal fluc- 
tuations of the order parameter A = |A| exp(iip). Such 
fluctuations result in a temporary local destruction of | A| 
accompanied by the phase slippage in the corresponding 
points of the wire. According to the Josephson relation 
V = ip/2e this process must cause a non-zero voltage 
drop across the superconducting sample thus bringing it 
into a resistive state. 

Quantitative theory of these thermally activated 
phase slips (TAPS) was worked out by Langer and 
Ambegaokar— and by McCumber and Halperin 4 . This 
LAMH theory predicts that in a superconducting wire 
TAPS are created with the rates r± which are defined 
by the activation dependence 



r± = B±e 



-SF ± /T 



(1) 



Here 5F± are effective free energy barriers which the sys- 
tem should overcome in order to create a phase slip cor- 
responding to the overal phase change ±2ir. These po- 
tential barriers are essentially controlled by the super- 
conducting condensation energy for the volume of the 
TAPS core where the order parameter |A| gets destroyed 
by thermal fluctuations. In the absence of any external 
bias one naturally has T + = T_ and, hence, no net volt- 
age across the sample can occur. Applying an external 
current / one lifts the symmetry between "positive" and 
"negative" TAPS. As a result, there appears a voltage 
drop (and, hence, non-zero resistance R(T)) proportional 
to the difference between the two TAPS rates T + — r_. 

According to Eq. TAPS remain significant only at 
temperatures close to Tc while R(T) decreases exponen- 
tially as T is lowered well below the critical temperature. 



This prediction was fully confirmed in experiments^ 
where the activation behavior of R(T) was detected in 
small superconducting whiskers with diameters ~ 0.5 
fim. Later is was realized that in thinner wires not 
only thermal but also quantum fluctuations of the or- 
der parameter (quantum phase slips) become impor- 
tant. It was demonstrated both theoretically 7 -^ and 
experimentall y 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 that quantum phase slip 
effects can yield appreciable resistivity of superconduct- 
ing wires with diameters in the range ~ 10 nm even well 
below Tc- For more details we refer the reader to the 
review- 
Turning again to thermal fluctuations near Tc, we note 
that, while for evaluation of the free energy barriers for 
TAPS 6F± in Eq. (JTJ it suffices just to solve the standard 
Ginzburg-Landau (GL) equations, the problem of finding 
the pre-exponent B± is in general much more involved, 
as it requires employing the formalism which properly 
accounts for dynamical effects in superconductors. Mc- 
Cumber and Halperin (MH) 4 treated this problem within 
the formalism of the so-called time dependent Ginzburg- 
Landau (TDGL) equations 1 ^ which was available at that 
time. Unfortunately, this formalism is known to suffer 
from serious drawbacks (see, e.g., Refs. [l1 fl~7lll8l for fur- 
ther discussion) and it is in general hardly applicable 
below Tc- Thus, although MH calculation 4 - of the pre- 
exponent B± was correct and sound by itself, their final 
result needs to be re-analyzed on the basis of a more solid 
theoretical approach. This task will be accomplished be- 
low. 

The structure of the paper is as follows. In Sec. 2 
we briefly re-capitulate the microscopic effective action 



formalisi 



and employ it in order to estimate the 



fluctuation correction to the order parameter of ultra- 
thin superconducting nanowires. In Sec. 3 we apply this 
formalism in order to evaluate the pre-exponent B± in 
the expression for the TAPS rates (Q}. We will then per- 
form a detailed comparison between our result and that 
of Ref. 4. In addition, we will present a simple for- 
mula which expresses the voltage noise in superconduct- 
ing nanowires via the TAPS rates ([I]) evaluated here. 
Some technical details of our calculation of fluctuation 
determinants are relegated to Appendix. 
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II. EFFECTIVE ACTION AND GAUSSIAN 
FLUCTUATIONS 

Consider a uniform superconducting wire with cross 
section s and length X. In order to account for super- 
conducting fluctuations in such a wire we will use the 
effective action approach developed in Refs. HUll Our 
starting point is the path integral representation of the 
grand partition function 



VAVVVA e~ Srff , 



(2) 



where S c g [A, V, A] is the imaginary time version of the 
effective action for a superconducting wire. The fluctu- 
ating order parameter field A as well as the scalar and 
vector potentials V and A depend on coordinate x along 
the wire (i.e. —X/2 < x < X/2) and imaginary time 
r restricted to the interval < r < 1/T. The exact 
expression for this effective action is obtained by inte- 
grating out the electron degrees of freedom and is not 
easily tractable in a general situation. In order to sim- 
plify this general expression for the action S e g one can 
assume that deviations of the amplitude of the order pa- 
rameter field A(x, r) from its equilibrium value Ao are 
relatively small. This assumption allows to expand the 
effective action in powers of SA(x, t) = A(x, t) — Ao 
and in the electromagnetic fields up to the second order 
terms. The next step is to average over disorder. After 
such averaging the effective action becomes translation- 
ally invariant both in space and in time. Performing the 
Fourier transformation we obtain&ii 

diodq f\Af_ + C\V\ 2 
Ls s 



<7 - £ 



Xd 



qV 



-A 

c 



XJ 



iuj 



XL 

Am 2 



iqtp- 



—A 

c 



Xa\5A\ 



(3) 

The functions XA^,q), Xj(v>Q), Xl(^,<?) and xr>(w,<?) 
are related to the averaged products of the Matsubara 
Green functions. The corresponding general expressions 
arc established in Refs. [HQS They are rather cumber- 
some and will not be specified here. In what follows we 
will use only simplified forms of these functions applica- 
ble in certain limits. 

As the action S e s © is quadratic both in the voltage 
V and the vector potential A, these variable can be inte- 
grated out exactly. After such integration one arrives at 
the effective action S which only depends on ip and 5 A. 
We get 



dtudq 



{H^,<1)W\ 2 + Xa\5A\ 2 } 



(4) 



Since usually the wire geometric inductance L remains 
unimportant, in what follows we will disregard this quan- 
tity by setting L — 0. Then we obtain^ 



' XJ , ,2 
>4? W 



XL 

4:171* 



I 2 ) C- 



xoq 



X.iXl „2 



+ x.j + XDq 2 



(5) 



The effective action ((4]) allows to directly evaluate the 
fluctuation correction to the equilibrium value of the or- 
der parameter in superconducting nanowires. Performing 
Gaussian integration over both <p and SA we arrive at the 
wire free energy 



F = F> 



bcs — tt 



T \ - 

2 



in 



^F{uJ,q) Axa(w,<?) 
2^V A 2 + 11 ' 



2N 



(0) 



where Fbcs is the standard BCS free energy and A is the 
BCS coupling constant. The order parameter is defined 
by the saddle point equation dF/dA = and can be writ- 
ten in the form A = Ao — SAq , where Ao is the solution 
of the BCS self-consistency equation OFbcs /dAo = 
and the fluctuation correction SAq has the form 



SA = - 



T f d 2 F BCS Y 



OA 2 



dA, 



E 



In 



J 

A-T 7 ^, q) 
2N A 2 



hi 



AxA(^,g) 



2A^o 



(7) 



First let us consider the low temperature limit T <C 
Ao . It is useful t o note that at large values of the wave 
number \q\ ^ Aq/D and/or frequency \u>\ 3> Ao the 
functions T{uj, q)/A 2 and^A^,?) are weakly affected by 
superconductivity. Hence, we can restrict the sum in Eq. 
(J7|) only to low frequencies \w\ < A and wave numbers 
\q\ < yj 'Aq/D. It will be convenient for us to introduce 
dimensionless parameters y = u/Aq, z = q*J D / Aq and 
express the kernels as follows 

XA = N F A (y,z); X J = e 2 N Fj(y, z); 

e 2 N n D 



Xl = m N a DA F L (y,z); xd = 



-F D (y,z);(8) 



where all the functions Fj are dimensionless. The func- 
tion T^^q) acquires the form 



N A 2 



(y 2 Fj + z 2 F L )(C/ S + e 2 N Q z 2 F D ) 
C/s + e 2 N Fj + e 2 N z 2 F D 
e 2 N z 2 FjF L 



C/s + e 2 N Fj + e 2 N z 2 F D 



(9) 



For a wire of length X we obtain dF/dAo = 2NqsX and 
at T — the correction to A reads 



<5A 



The integral 



A 



8sN y/DAo 
Fa(v,z) 



dydz 



Ill 



(10) 



dydz F A (y,z) 



In 



(2tt)2 F A (1,1) 

is well convergent at small y and z, therefore we can 
replace it by a constant of order one. The integral 

dydz F{y,z) 
(2tt)2 V(l,l) 
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is only slightly more complicated, since T(y, z) — > for 
y, z — » 0. However, since the function T(y, z) enters only 
under the logarithm, this integral is convergent as well. 
Making use of the above expressions for the functions 
T[u), q) and xa(w, q) at T — > we obtain 



A 



1 

9i 



Gi 



3/2 
ID ■ 



(11) 



Here is the dimensionless conductance of the wire 
segment of length £ and Giro is the Ginzburg num- 
ber for a superconducting nanowire defined as the value 
(T<7 — T) /Tc at which the fluctuation correction to the 
wire specific heat becomes equal to the specific heat jump 
at the phase transition point. In the case of quasi- ID 
wires this number readsi^: 



Gi 



0.15 



ID 



(12) 



We note that in Eq. (TTTj) fluctuations of both the phase 
and the absolute value of the order parameter give con- 
tributions of the same order. The estimate (|lip demon- 
strates that at low temperatures suppression of the order 
parameter in superconducting nanowires due to Gaus- 
sian fluctuations remains weak as long as ^ > 1 and 
it becomes important only for extremely thin wires with 
GiiD ~ 1 in which case the width of the fluctuation re- 
gion ST is comparable to Tc and the BCS mean field 
approach becomes obsolete down to T = 0. 

Turning to higher temperatures we observe that at T 
sufficiently close to the critical temperature Tc it is nec- 
essary to retain only the contribution from zero Matsub- 
ara frequency. At the same time the terms originating 
from all non-zero frequencies are small in the parame- 
ter Ao (T)/Tc <C 1 and, hence, can be safely omitted. 
Performing the integration over q we get 



^Aq 
Ao(T) 



T 

5F' 



where 



Ao(T) = 



l8ir 2 T(T c -T) 



7C(3) 



and 



21C(3)' 



T) 



3/2 



(13) 



(14) 



(15) 



turns out to be exactly equal to the magnitude of the 
effective free energy barrier for TAPS in the LAMH 
theory 3 "^ in limit of small transport currents (see below). 
Eqs. (TIB"]) , (fl"5|) demonstrate that at temperatures close 
to Tc Gaussian fluctuations of the superconducting or- 
der parameter in thin wires become more significant and 
effectively wipe out superconductivity at SF < Tc, i.e. 
already in much thicker wires than in the case of low 
temperatures T -C Tc- 



III. THERMALLY ACTIVATED PHASE SLIPS 

In what follows let us restrict our attention to super- 
conducting wires in which the condition SF 3> Tc is well 
satisfied and, hence, the effect of Gaussian fluctuations 
on the order parameter A (T) can be safely neglected. 
This condition requires the wire to be sufficiently thick 
and/or the temperature should not be too close to Tc, i.e. 
(Tc — T)/Tc 3> Gim- At the same time we assume that 
the temperature is still not far from Tc, i.e.Tc — T<^ Tc 
in which case the physics is dominated by thermally acti- 
vated phase slips 3 -^. As we already discussed, sufficiently 
thin superconducting wires acquire non-zero resistance 
even below Tc due to TAPS, and this resistance is essen- 
tially determined by the TAPS rates ((TJ. 



A. Activation exponent 

The free energy barriers SF± for TAPS corresponding 
to overall phase jumps by ±2tt have been evaluated by 
Langer and Ambegaokar—. Here we briefly re-capitulate 
their results. In order to obtain 6F± entering into Eq. 
(fTj) we make use of the standard Ginzburg-Landau free 
energy functional for a wire of length X: 



F[A(a 



= siVo / ax 

-X/2 



ST 



dA 



dx 



T-T c 
T C 



7C(3) , 
16tt 2 T 2 1 



2e 



[cp(X/2) - <p(-X/2)]. 



(16) 



Here ip(x) is the phase of the order parameter A(x) and 
I is the external current applied to the wire. 

The saddle point paths for this functional are deter- 
mined by the standard GL equation 



nD d 2 A T-T, 



c 



8T dx 2 



T c 



A 



7C(3) 
8vr 2 T 2 



A| 2 A = 0. 



For any given value of the bias current 



ttbNqDs , A |2 



2T 



A^V^ 



(17) 



(18) 



this equation has a number of solutions. The TAPS free 
energy barrier 5F+ is determined by the two of them. 
The first one, A m = |A m | exp(i<^> m ), corresponds to a 
metastable minimum of the free energy functional. This 
solution reads 



,1 + 2 cos a 

Ao(T)W , (p rn 



2T 



Ix 



neN Ds |A m | 2 



.(19) 



Here Ao(T) is the equilibrium superconducting gap de- 
fined in Eq. (|14p and the parameter 



3 Ic (1 - 2(/// c ) 2 ) V ' 
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accounts for the external bias current I. The Ginzburg- 
Landau critical current Iq is defined by the standard 
expression 



The second, saddle point, solution A s (x) 
| A s | exp(i(p s ) of Eq. (fT7|) has the form 



(21) 



Ao(T) 



\ 



1 + 2 cos a 



2 cos a — 1 



2T7 



dx' 



cosh 
dx' 



V2cosa — 



A iWo |A s (x')| 2 ' 



(22) 



where £(T) = ^JirD/A(T c — T) is the superconducting 
coherence length in the vicinity of Tc- 

The free energy barrier SF + 2tt in Eq. |T]) is set by the 
difference 



SF+ =F[A s (x)} - F[A m (x)] = SF 



V2 cos a — 1 



2 I ( V3 /2cosa - 1 

arctan — i / 

3J C V 2 V 1-cosa 



(23) 



where SF is defined in Eq. (|15p . The free energy barrier 
for "negative" TAPS is determined analogously and 
is related to SF + as follows 



irl 

SF- = SF+ + — . 

e 



B. Pre-exponent 



(24) 



Now let us turn to the pre-exponent B±2 7r in the ex- 
pression for the TAPS rate ([T]). For simplicity we first 
analyze the TAPS rate in the zero current limit in which 
case SF + = <5F_ = SF and B + = B- = B. In order 
to evaluate B one should go beyond the stationary free 
energy functional p^|) and include time-dependent fluc- 
tuations of the order parameter field A(x,t). In Ref. 4 
this task was accomplished within the framework of a 
TDGL-based analysis. Employing TDGL equation it is 
possible to re-formulate the problem in terms of the cor- 
responding Fokkcr-Planck equation 1 - which can be con- 
veniently solved for the problem in question. Since the 
important time scale within the TDGL approach is the 
Ginzburg-Landau time 



T~GL 



i\T c -T\ 



(25) 



this time also naturally enters the expression for the pre- 
exponent B derived in Ref. 0. 

Unfortunately the TDGL approach fails below Tc- For 
the sake of illustration, let us for a moment ignore both 



the scalar and the vector potentials. The TDGL action 
for the wire is then usually written in the form 



Stdgl = N Ts 



E 



dx^|A| 2 
8T 1 1 



c2tF[A(x,t)L26) 



where the GL free energy functional F[A(x,t)] is de- 
fined in Eq. (|16[) . This form can be obtained from the 
action ^ by formally expanding the kernel xa in Mat- 
subara frequencies and wave vectors uj^,Dq 2 -C AwT. 
Note, however, that since the validity of the GL expan- 
sion is restricted to temperatures T ~ Tq, the Matsub- 
ara frequencies \lu„\ — 2n\n\T cannot be much smaller 
than 4irT for any non-zero n. Hence, the expansion 
#(1/2+ |w„|/47rT) - #(1/2) -> 7r|w„|/8T - which yields 
TDGL action (f2"o| - is never correct except in the sta- 
tionary case u> n = 0. Already these simple arguments 
illustrate the failure of the TDGL action (|26|) in the Mat- 
subara technique. Further problems with this TDGL ap- 
proach arise in the presence of the electromagnetic po- 
tentials V and A. We refer the reader to the paper— for 
the corresponding discussion. 

In view of this problem one should employ a more ac- 
curate effective action analysis. Since the microscopic 
effective action for superconducting wiresi&il cannot be 
easily reduced to any Fokker-Planck-type of equation 
it appears difficult to directly employ the McCumber- 
Halperin approach^ in order to evaluate the pre-exponent 
B in the expression for the TAPS rate ([T]). For this rea- 
son, below we will proceed differently and combine our 
effective action formalism with the well known general 
formula for the decay rate of a metastable state expressed 
via the imaginary part of the free energy. This method 
is applicable provided the system is not driven far from 
equilibrium. For the decay rate in the thermal activation 
regime one ha a 20 ' 21 i 22 



r = -2— ImF(T), 



(27) 



where T* is an effective crossover temperature between 
the activation regime and that of quantum tunneling 
under the potential barrier. Formally T* is defined as 
temperature at which a non-trivial saddle-point solution 
A(t, x) describing QPS first appears upon lowering T. 
Within the accuracy of our calculation it is sufficient to 
estimate T* simply by setting the QPS action Sqps(T) 
equal to activation exponent, i.e. 



S QPS (T*) ~ SF{T*)/T* 
At sufficiently small currents one has^ 



Sqps{T) = AsN oy /N A {T), 



(28) 



(29) 



where A is a numerical constant of order onei Hence, the 
condition (j55J) yields A (T*) - T* or, equivalently, T* = 
aTc, where the numerical factor a < 1 is sufficiently close 
to unity, i.e. T* ~ T c . As the whole concept of TAPS 
is only valid at T close to Tc, one always has T* /T ~ 1. 
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Thus, with the same accuracy one can actually use the 
expression for the decay rate in the quantum regim e 22 ! 23 
r = -2ImF(T), c f. Ref. Q. Here we will retain the 
parameter T* for the reasons which will be clear below. 

Following the standard procedure we expand the gen- 
eral expression for the effective action around both so- 
lutions (fT5|) and (|2"2")l up to quadratic terms in both the 
phase tp and the amplitude SA. One can verify that in 
the limit Aq(T) <C T the contributions from fluctuating 
electromagnetic fields can be ignored and we obtain 



S s/m =F{A s/m ]+S 2 S. 



s/m-i 



(30) 



where 



x [5A(uj n ,x)xA l \\u n \;x - x')dA(u n ,x') 
+ Lp(uj n ,x)k { ™ l) (\uj n \\x - x')ip(uj n ,x')], 

S 2 S S = — J dxdx 1 

x [&A(uj ni k)xa (knh x, x')SA(uj n , x') 

+ ip(w n ,x)kg>(\u n \;x,x')cp(GJ n ,x')]. (31) 

Here u n = 2-kTu are Bose Matsubara frequencies. The 
functions Xa"^ an d are expressed in terms of the 

kernels xa, XJ an d XL as follows: 

X { ?H\u; n \;x-x')= f 



2vr 



kW(\u> n \;x-x')= I g e ^-') 



The functions 



Is) 

and kin describe fluctuations around 



the coordinate dependent saddle point A s (x), and, there- 
fore, cannot be easily related to xa, XJ an d XL- Fortu- 
nately, the explicit form of Xa an d 4 1Sl not important 
for us here. 

The pre-exponent B in eq. |T]) is obtained by integrat- 
ing over fluctuations 8 A in the expression for the grand 
partition function. One arrives at a formally diverging 
expression which signals decay of a metastable state. Af- 
ter a proper analytic continuation one finds the decay 
rate in the form JT]) with 



B 



-2T* Im JJ 



detxA^n) detfc^ m) (u;n) 
detxA^n) detfc^ s) (w„) 



(32) 



Here it is necessary to take an imaginary part since one 
of the eigenvalues of the operator (0) is negative. 

The key point is to observe that at T ~ T c all Matsub- 
ara frequencies \u> n \ = 2i:T\n\ - except for one with n = 



- strongly exceed the order parameter, \ui n \ ^> A (T). 
Hence, for all such values the function XA^ml) ap- 
proaches the asymptotic form&i£ which is not sensitive to 
superconductivity at all at such values of u) n . Hence, as 
long as Aq(T) «Twe have detxA^( w n) 



and detkjp\cj n ) ~ det k^' 1 ' (ui n ). The corresponding de- 



detxA l) (u n ) 



terminants in eq. (|32|) cancel out and only the contribu- 
tion from uj n — remains. It yields 



B 



-2T* Im 



detx^W det4 m) (0) 
\\ detxi s) (0) det4 s) (0) ' 



(33) 



The ratio of these determinants can be evaluated at zero 
current with the aid of the GL free energy functional (JTJ 
This calculation yields, see Appendix A: 



Im , 



detxA^O) det4 m) (0) 
\ detx^W det4 s) (0) 



M * /E (34) 



where, as before, SF is defined in Eq. (fT5)) . 

Combining the above expressions we arrive at the final 
result for the TAPS rate in the zero bias limit: 



r _ r 4^ X [SF 
Z(T) V T 



5F_ 

T 



(35) 



Turning to the case of non-zero bias one can essentially 
repeat the whole calculation which now yields two differ- 
ent TAPS rates T±. Of practical importance is the limit 
of transport currents I sufficiently close to the critical 
one, i.e. 1 — I/Ic *C 1- In this regime T_ is negligibly 
small whereas T + , on the contrary, increases since the 
free energy barrier 



SF+(I) 



4 x 6 3 / 4 
15" 



-SF 1 



i) 



5/4 



(36) 



becomes lower than that at smaller currents. Accord- 
ingly, TAPS can be detected easier in this limit 2 ^. 

The pre-exponent B + has essentially the same form 
as that defined by Eq. (f3"5]) . one just needs to replace 
SF -> SF+(I) and T* -> T*(I). In the limit T c -T < T c 
considered here the current dependence of the crossover 
temperature T* appears insignificant in most cases and 
with sufficient accuracy one can set T*(I) ~ T*. Indeed, 
very generally one can express T*(I) = T* f(I/I c (T*)), 
where Ic{T*) is the critical current at temperature T* 
and f(x) is some universal function with f(x -C 1) — 1. 
Having in mind the strong temperature dependence of 
Ic(T) in the temperature interval Tc — T <C Tb we find 

I/I C (T*) < I C (T)/I C (T*) ~ {T c - Tf/ 2 /T c /2 « 1, 

and, hence, T*(I) » T*(0) = T*. Thus, in the vicinity 
of the critical current I C (T) - 1 < I C (T) the TAPS rate 
can be expressed in the form 



8.84T* 
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where SF + (I) is defined in Eq. (|36[) and the numerical 
prefactor is again established from the calculation of the 
fluctuation determinants which is fully analogous to that 
presented in Appendix A. 

Summarizing all the above results and substituting 
T* = aTc we arrive at the final expression for the TAPS 
rates 



r ± (/) = KaT c ^/^)exp 



T 



SF ± (I) 
T 



(38) 



where k(I) is a smooth function of I varying from k(0) ~ 
5.53 to k — 8.74 at I c — I <C Ic an d, as before, the 
numerical prefactor a is of order (and slightly smaller 
than) one. Eq. (|38[) is the central result of this work. 
This expression is supposed to be valid at Tc — T <C Tc 
and at any bias current / < Ic as long as 5F±(I) ^> T . 



C. Comparison with McCumber-Halperin result 



I. This analysis yields 



T, 



TDGL 



TDGL 



I 

1 

Ic 



5/4 



(41) 



Combining Eqs. ((40 
at the pre-exponent 



Btdgl 



4T1) with the result (l37l) , we arrive 



15/8 




TGL C(T) 



which is again in the agreement with Ref. 4. Thus, with 
the aid of the general formula (f2"T)) describing thermally 
activation decay of a metastable state we confirm that 
the McCumber-Halperin result for the TAPS rate is es- 
sentially correct within the TDGL-type of formalism. Un- 
fortunately, however, the latter formalism is inaccurate 
by itself. In particular, in the expression for the TAPS 
rate it does not allow to correctly obtain the classical-to- 
quantum crossover temperature T*. 



Let us compare our result (|38l) with the expression for 
the TAPS rates derived in Refill from the TDGL-type of 
analysis. We observe that Eq. (|38|) does not contain the 
Ginzburg-Landau time tql and exceeds the correspond- 
ing expression!, by the factor ~ T*t gl ~ {l~T/T c y l > 
1. On top of that, in the vicinity of the critical cur- 
rent the pre-exponent in the TAPS rate (|37|) depends 
on I as B + oc (1 — I/Ic) 5 ^ s in contrast to the result^ 
Btdgl oc (1 - I /I C ) 15/8 - 

In order to understand the origin of these differences let 
us - just for the sake of illustration - for a moment adopt 
the TDGL action {26]) and re-calculate the TAPS rate 
Ftdgl employing Eq. (|27|) . Since the whole calculation 
of the fluctuation determinants remains the same (see 
Appendix A) we should only re-evaluate the crossover 
temperature which we now denote as T^ DGL . To this 
end we again first set I — > and consider fluctuations of 
the order parameter around the saddle point A s (x) along 
the unstable direction (|A15p choosing 



5A(t, x) = i 



cos (2itT£ dgl t) 
\72cosh(a;/£) 



a 



(39) 



where C is a constant. Substituting this expression into 
the linearized TDGL equation and formally treating tql 
as an independent parameter, we define the classical-to- 
quantum crossover temperature T^ DGL as that at which 
a non-zero solution (C ^ 0) first appears. This definition 
yields 



T TDGL = 



(40) 



D. Temperature-dependent resistance and noise 

In order to complete our analysis let us briefly address 
the relation between the above TAPS rate and physi- 
cal observables, such as, e.g., wire resistance and voltage 
noise. Every phase slip event implies changing of the su- 
perconducting phase in time in such a way that the total 
phase difference values along the wire before and after 
this event differ by ±27r. Since the average voltage is 
linked to the time derivative of the phase by means of 
the Josephson relation, (V) = (<p/2e), for the net voltage 
drop across the wire we obtain 



V = -[T + (I)-T.(I)] 

e 



(42) 



where T± are given by Eq. J38|) . In the absence of any 
bias current I — > both rates are equal T± = T and the 
net voltage drop V vanishes. In the presence of small 
bias current I <ti Ic we obtain 



r±(/) = r e 



_ r ±Kl/2eT 



(43) 



Thus, at such values of I and at temperatures slightly 
below Tc the I — V curve for quasi- ID superconducting 
wires takes a relatively simple form 



V 



— r sinh , 

e 2eT' 



(44) 



The zero bias resistance R(T) = (dV/dI)i—o demon- 
strates exponential dependence on temperature and the 
wire cross section 



Substituting (|40|) into Eq. (f35|) we arrive at the expres- 
sion for Ttdgl just 2 times bigger than that derived in 
Ref. in the limit / ->■ 0. 

An analogous - though slightly more complicated - 
analysis can be performed also at non-zero bias current 



l R{T) 
2tt 



2V6tt 



aT c X 
T £(T) 




exp 



SF 



(45) 



To complete our description of thermal fluctuations 
in superconducting wires we point out that in addition 
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to non-zero resistance (|45|) TAPS also cause the voltage 
noise below Tq- Treating TAPS as independent events 
one immediately concludes that they should obey Pois- 
sonian statistics. Hence, the voltage noise power 

S v = 2 / dt(SV(t)SV(0)) 



is given by the sum of the contributions of both "positive" 
and "negative" TAPS, i.e. 



s v = ^-[r + (i)+r.(i)] 



(46) 



At small currents I Ic this expression reduces to the 
following simple form 



47r 2 irl 
Sv = — Tcosh— . 



(47) 



Similarly to the wire resistance the voltage noise rapidly 
decreases as one lowers the temperature away from Tq- 
Only in the vicinity of the critical temperature this TAPS 
noise remains appreciable and can be detected in exper- 
iments. 

In conclusion, we have demonstrated that the rate 
for thermally activated phase slips in superconducting 
nanowires evaluated within the microscopic effective ac- 
tion analysis turns out to be parametrically bigger as 
compared to the TAPS rate derived from the TDGL-type 
of approach. Simultaneous measurements of both TAPS- 
induced resistance and noise appears to be an efficient 
way for quantitative experimental analysis of thermally 
activated phase slips in superconducting nanowires. 
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APPENDIX A: EVALUATION OF 
FLUCTUATION DETERMINANTS 

Let us set 1 = and write the Ginzburg-Landau free 
energy in the form 



F = 



3SF 



X/2£ 



dr) 



2 2 
— U —V 



{u 2 + v 2 f 



IA1) 



-x/n 



where we introduced r\ = x/£, u = ReA/Ao(T), v = 
ImA/A (T). In terms of these dimcnsionlcss variables 
the metastable solution (1191) reads 



u m (r)) = 1, v m {rj) = 0, (A2) 
while the saddle point solution takes the form 

u a {ri) = tanhr/, v s (rj) = 0. (A3) 



The second variation of the free energy (|Aip around 
its saddle point with v(y) = is 



X/2£ 



Su' 2 + Sv' 2 

" + r > - Su 2 - Sv 2 



-X/2£ 

+u 2 (Su 2 +5v 2 ) + 2u 2 Su 2 



(A4) 



Here Su and Sv describe fluctuations of respectively the 
absolute value and the phase of the order parameter. Ac- 
cordingly the operators \A and k v read 



xi s) (o) 



xi m) (o) 



^ s) (o) 

fcW(O) 



3SF 


' d 2 




6 


~4T 


dr] 2 


+ 4 


cosh 


35F 


' d 2 






~W 


dr/ 2 


+ 4 




35F 


- d 2 




2 


~4T 


dr\ 2 


cosh 2 r]_ 


35F 


' d 2 ' 






~4T 


dr] 2 







(A5) 



In order to fix the boundary conditions we note that fluc- 
tuations of the absolute value of the order parameter in 
the bulk leads are negligible. Hence, we set 



5u(-X/2) = Sv(X/2) = 



(A6) 



Likewise, since the current density vanishes in the bulk 
leads we can choose 



5v'{-X/2) = 8v'{X/2) = 0. 



(A7) 



Let us evaluate the eigenvalues of, say, the operator 
X ( a(0)- These eigenvalues Al s) = (35F /AT)X ( n s) are ob- 
tained from the Schrodinger equation 



dr] 2 



6 



cosh T] 



Su = (A - 4)6w 



(A8) 



with appropriate boundary conditions. The correspond- 
ing localized solutions of this equation have the well 
known form 2 ^: 



A 



W _ 



0, Su{ s) ( V ) = 



1 



4 cosh 2 T] 
3 sinh?7 



2 cosh 2 



n 



(A9) 



In order to find the eigenvalues in the continuous spec- 
trum we introduce transmission, £(A), and reflection, 
r(A), amplitudes of the potential well 2 ^ 



t(A) = 

V 1 (l + zVA^4)(2 + iVA^4) V ' 



0. (A10) 



Thus, at large negative r] the wave function has the 
form 8v(rf) = Cie 1 ^^^ + C 2 t{\)e- l ^ :zii ^ , while at 
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large positive r/ the same wave function is 8v(rf) = 
Cit(A)e iv/A ~ 4r ' + C2e~ 4%/A ~ 4 ' 7 . Imposing the boundary 
conditions (|A6|) we arrive at the following equation for 



the eigenvalues A 



(n = 3,4,. 



where 



^(A) = - U(A)e 



t(X) 



(All) 



(A12) 



In the limit X/£ — ► 00 Eq. (|A11[) also applies for the 

1. 



discrete eigenvalues A^, A^ 



The eigenvalues of the operator Xa (0) are obtained 
analogously. They are defined by the equation 

F m {\W) = 0, n=l,2,... (A13) 



where F m (X) = sin (yA — 4 A/£). Observing that 
i^s.mCA) oc n^Li(A~ A^ m ) and extracting the zero eigen- 
value of Xa (0) m a standard way, after the integration 
over this zero mode we obtain 



/ detxg(O) 
detxl(O) 



The 



ratio 



UX 
3? 

of 



3 SF 



lim lim 



4 2ttT a^ox- 





fluctuation 



determinants 



operator kp(0) has two localized eigenf unctions fa^fa) 
and <5^ s) (77) with the eigenvalues e[ s \ = (3SF/4T)e{% 



1 1 

cosh 77 ' 



(A15) 



e 2 s) =0, = y^tanhr;. (A16) 



(a) 

The eigenvalue is negative and, as usually, it is asso- 
ciated with the unstable direction in the functional space. 
The ratio of the determinants is expressed as follows 



/ detAff(O) _ / Cm (A) 
detfcj(O) aTo\/G s (A)' 



(A17) 



where G S (A) = sin \y/\X/£\, while G S (A) reads 



1 / -iV\4- 



(A18) 



where t(A) = (i\/A- !)/(! + «VA). Thus, we get 



/detAff(O) 
detfc-(O) 



(A19) 



detfc™(0)/detfcj(0) is evaluated analogously. The Combining Eqs. (|A14|) and (|A19|) . we arrive at Eg. ([M]) . 
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